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Abstract 

The set of integrable symiiietries of the nonstationary Schrodinger equation 
is shown to admit a natural decomposition into subsets of mutually commuting 
symmetries. Hierarchies of time evolutions associated with each of these sub- 
sets ultimately lead to nonlinear (possibly, operator) equations of the Kadomtsev- 
Petviashvili I type or its higher analogues, thus demonstrating that the linear prob- 
lem itself constructively determines the associated nonlinear integrable evolution 
equations and their hierarchies. 
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1 Introduction 

In this article, we use the nonstationary Schrodinger (NS) equation as an example to 
demonstrate that a given linear problem determines the nonlinear evolutions associated 
with it. Our approach is based on the study of the set of symmetries for the NS equation 
using the inverse scattering transform (1ST). We here consider those symmetries that can 
be called integrable, which means that the associated flows can be expressed in terms of 
the potential of the linear problem or in terms of some special values of the Jost solutions. 
The set of all integrable symmetries is decomposed into subsets of mutually commuting 
ones; with each of these subsets, we associate a hierarchy of evolution equations. We show 
that all these subsets of commuting symmetries have a natural hierarchy structure, and 
we investigate the corresponding evolutions. We find that the Jost solution of the original 
NS equation with respect to the two lowest times of any such hierarchy also obeys an NS 
equation with an operator-valued potential. Moreover, the third (or any higher) time of 
the corresponding hierarchy provides an operator differential equation whose linear part 
coincides with the linear part of the Kadomtsev-Petviashvili I (KPI) equation y (or of 
its higher analogue respectively). This linearized KPI equation naturally appears as the 
equation for the spectral data here; it follows from a simple commutator identity that 
the spectral data satisfy this equation (or its higher analogues) regardless of the choice 
of the subclass of commuting symmetries. We also show that there are two distinguished 
subclasses of commuting flows characterized by the property that the above-mentioned 
operator potential is just a multiplication operator. These subsets of symmetries corre- 
spond to the standard evolutions and the ones that can be considered dual to them. We 
consider these subsets in detail and show that the second one leads to a transformation 
of both the dependent and independent variables, which essentially extends the range of 
applicability of the 1ST. 

The study of different symmetries associated with the NS equation (or other linear 
problems) goes back to the dressing procedure |g, || and was investigated either by these 
methods and their extensions or in formal algebraic terms (see |4|]- [p|). The approach 
adopted below is peculiar: we work strictly in the 1ST framework |l(]|]-[|l3[, starting with 
the NS equation with respect to (xi, 2:2) variables with a real potential u{xi,X2) assumed 
to be a smooth function that rapidly decays at infinity and satisfies some weak conditions 
on its norm |14| to guaranty the solvability of the direct and inverse problems. 

Our approach is based on a special version of the 1ST, the so-called resolvent approach 
(see [|l3|-[|9|)- In the next section, we briefly describe some constructions that enable us 
to combine the methods of the 1ST with the algebraic schemes of the type in [||-0. To 
control the transformations of all the involved objects — spectral data, Jost solutions, and 
the potential — we use the so-called scattering theory on a nontrivial background |20) . In 
Sec. 3, we consider integrable symmetries in this framework and derive general results 
related to their commuting subsets (evolutions). The simplest examples of such evolutions 
are given in Sec. 4. 

2 Extended operators and resolvent approach 

The first step in the resolvent approach consists in a special extension of the differential 
operators. Precisely, any given differential operator A with the kernel 

A{x,x') = A{xi,X2,d^,,dx.2)S{xi - x[)S{x2 - x'2) (2.1) 



is replaced with the differential operator A{q) with the kernel 

A{x, x'- q) = e"«(---')A(a;, x') = A{x^,X2, d^, + qi,d^, + q2)6{x - x'), (2.2) 

where qx — qiXi + q2X2 (we consider only the two-dimensional case here). In ( ^.2| ), q is 
a real (two-component) vector; below, we often omit explicit indication of the operator 
dependence on the variable q. The standard definitions of the dual operator, the complex 
conjugate operator, and the Hermitian conjugate operator are modified as follows: 



^'^"'''(a;, x'; q) = A\x' , x; -q), A*{x, x'; q) ^ A{x, x'] q), 

AHx,x';q) = A'ix',x;-q), (2.3) 

where the superscript t denotes transposition in the matrix case. In what follows, we 
refer to the kernel A(x, x'; q) of the operator A as the kernel in the x representation to 
distinguish it from the kernel of the same operator in the (p, q) representation defined 
via the "shifted" Fourier transform. 



A{p;q) = j^ J dx J dx'e'^''+^'^^^-'^-^' A{x,x';qc,), 



(2.4) 



where p and q are respectively real and complex two-dimensional vectors. In these terms, 
the meaning of this extension of the kernels of differential operators can be clarified by 
observing that the kernels A{p; q) of differential operators are polynomial functions of 
the variable q. As shown in iQ, this analyticity property of the extended differential 
operators considerably simplifies the study of their spectral transform. 

Below, we consider integral operators A{q) {A, for short) with kernels A{x,x';q) be- 
longing to the space S' of the Schwartz distributions with respect to all six real variables. 



xi, X2, x'l, x'2, qi, and (72 (or in the space of the {p, q) variables with respect to (2^)). We 
also consider a special subclass M C S' oi such operators, whose kernels are Schwartz 
distributions with respect to the variables x and x' and are piecewise continuous functions 
with respect to the variable q. 

The relevance of this subclass is demonstrated by the following simple example. Let 
Dj denote the extension of the differential operator dxj , j = 1, 2, by (^.2|), i.e., 

D,ix,x';q) = id^^+q,)Six-x'), J = 1,2. (2.5) 

This operator is uniquely invertible on the subclass A4 . Its inverse is given by the kernel 

D-\x,x';q) = sgnq,e-'^^^^^-^'^'^e{qj{x^ - x'^))5{x.j+^ - a^^+j), j = 1,2, (2.6) 

{i + 1 is understood mod 2) and is exactly the standard resolvent of the operator dx . 
This observation gave the name resolvent approach for the method based on the above 



extension of the original differential operators. The uniqueness of inverse (2^) fol- 
lows because the only (up to a factor) annulator of Dj belonging to S' has the kernel 
S{q,)S{x^_^i — x'j_^_i), which does not belong to A^. We also mention that if the q depen- 
dence of the extended differential operator can be trivially removed by multiplying the 



kernel by expq{x — x') (cf. (2.2)), the q dependence of the inverse operator cannot be 



removed in this way, as is shown by the example in (2.(;). 
We turn to the NS operator 



Cx^idx,+dl~u{xi,X2). (2.7) 



Using ( |2.2| ), we can write its extension as 



f/, 



where Lq is the extension corresponding to the zero-potential case and is equal to 



Lq = iD2 



Di 



by (2^). In (p.q), we introduce the multiplication operator U with the kernel 

C/(x, x'; q) — u{x)S{x — x'). 



The main object in our approach is the (extended) resolvent of operator (2 
as the inverse of the operator L on the subclass A4, i.e., 



where / is the unity operator, 



LM = ML = I, 



I{x, x' \ q) = S{x — x'). 



(2.9) 

(2.10) 
defined 

(2.11) 



We mention that the kernel of the resolvent Mq of the zero-potential operator Lq 
is equal to 



(2.12) 
9|) 



Moix,x':q) = 



sgn(a:2 - X2) 

2TTi 



daO 



((92 - 2agi)(x2 - 4))e-[«+*^("+»9i)l(--')^ (2.13) 



where we introduce the special two- vector 



(2.14) 



Therefore, this resolvent indeed belongs to M., and it is easy to see that it is unique in this 
subclass. In general, the existence of an inverse in the space of Schwartz distributions with 
respect to the variables x and x' follows from the Hormander theory, but the necessary 
properties with respect to the variable g, as well as the uniqueness statement, must be 
proved. This goes beyond the frame of this paper. We therefore assume here that the 
potential u{x) is a smooth function rapidly decaying at infinity; the resolvent therefore 
exists and is unique in M.. The class of such potentials is not empty because imposing 
the same small norm assumption used in |14| is sufficient to prove the existence of the 
Jost solutions. Moreover, we know (see |l6|] ) that the resolvent in fact exists in much 
more complicated situations, for instance, when the potential has a nontrivial behavior 
at infinity. 

To demonstrate that the extended resolvent generalizes all notions of the standard 
approach, we mention that it is directly related to the Green's function of operator (2.7). 
Indeed, if we introduce 



G{x, x', k) = e("-"')^=('')M(a;, x'; £cj(k)), 
then it is easy to verify that 

L^G{x, x' ; k) = 5(x — x'). 



(2.15) 



(2.16) 



We already mentioned that L{p\ q) is an analytic (polynomial) function of q. Tak- 
ing into account that the standard product of any two operators in terms of the (p, q) 
variables takes the form (cf. (2.4)) 



{AB){p:q) = J dp'Aip-p';q + p')Bip';q), 



(2.17) 



we see that M{p; q) must in some sense be close to an analytic function of the variable 
q. As shown in psl Qq| , studying its lack of analyticity leads to the notion of the Jost 
solution. To be more precise, we introduce the operator d2Mo whose kernel in the (p, q) 
representation is equal to dMo{p; q)/9q2. Then it is easy to see that 



{d2Mo){x,x';q) 



1 



87r|gi 



exp< —ill 



I 92 



+ iq-^ ) {x - x') 



(2.18) 



\v){x, x'-q) = 26{x^ - x'^) j dq2(^{ML„)d2Mo) (x, x'; q), 
{uj\{x,x'-q)^25{x2-x'^) j dq^({d^M^){L^M)){x,x'-q). 



It was shown in jl^] that the lack of analyticity of the resolvent M is given in terms of 
two special reductions of the (truncated) resolvent itself. These operators are defined as 



(2.19) 
(2.20) 



That they are reduced values is more explicit in the (p, q) representation: for example, 
|i/)(p, q) = (MLo)(P7q)|q2=q2- Both kernels |i/)(p, q) and (w|(p, q) are independent of 
(\2 Slid are analytic functions of qi in the upper and lower half-planes with a possible 
discontinuity at the real qi axis. Both tend to 6{p) as qi -^ oo. We call these operators 
the Jost solutions of the NS equation (and its dual), because if we introduce the functions 



(2.21) 
(2.22) 

(2.23) 



X(a;,k) - / dxV'''*(^^-"i)|i.)(a;,x';kQ), 

C(a;',k) = /"da;e*'^''(^i-^i)(a;|(x,x';kcj) 

for some complex k, then the standard Jost solutions are given by 

$(x, k) = e-^^f'^^^xCa;; k), *(x, k) = e'^^^'^'' £,{x; k). 



We emphasize that $ and ^ are the Jost solutions of the original operator (pj), not of 
the extended one, i.e., they satisfy 



/:,<i>(x,k) = o, £^"-'vi/(a,,k) = o 



(2.24) 



while the e xtens ion o f the operator L causes the Jost solutions to be determined by 
reductions ( 2.1£ ) and ( 2.20 ) as functions of the complex spectral parameter k. Precisely, 



^d^^xi^, k) + dl^x{x, k) = u{x)x[x, k) + 2ik.d^^x{x, k). 



(2.25) 



and it is easy t o veri fy tha t by t he above definitions, x{^^ k) is normalized to 1 at k — > (X). 
Operators (2.19) and (2.2C) satisfy the "differential equations" 



L\v)^\u)Lo, (w|L = Lo(w|, 



(2.26) 



which are of course equivalent to Eqs. (2.24) after transformations (2.21), ( ^.22 ), and 
(^.25). The reality condition for the potential u{x) is equivalent to 



\t _ 



(^1, 



(2.27) 



which can also be written as $(x, k) = ^(x, k). 



In Pql, we also proved that the operators \v) and (ijj| are the inverses of each other, 

{u;\y)^I, \i^){uj\^I, (2.28) 

and that we have a bihnear representation of the operator L and of the resolvent M by 
their means, 



L = \v)La{uj\, 
M ^ \v)Mo{lj\ 



(2.29) 
(2.30) 



Using the resolvent approach, we thus obtain equations that have the explicit meaning 
of dressing the zero-potential operator Lq with the operators \v) and (wj. One feature of 
our a pproach is that these dressing operators are in turn given as reductions ( 2.19| ) and 
( 2.2C ) of the resolvent itself. 

To formulate the inverse problem, we use the analyticity properties of the kernels 
|j^)(p, q) and (a;|(p, q) and introduce their boundary values on the real axis from the 
upper and lower half-planes. In terms of the [x, q) representation, the corresponding 
kernels are defined by means of the limits 



\v^)(x,x';q)= ]iu\ \v){x,x'\q), 

{uj^\{x,x\q) = Vnn{uj\{x,x'\q), 
gi-+±0 



(2.31) 
(2.32) 



where the l.h.s.'s, of course, do not depend on q. Because the potential is real, they are 
mutually conjugate, 

|z.±)t = (c.T| (2.33) 

(see (2^) and ( 2.27 )). The spectral data are defined [^ Q as 



and by ( ^^ 



F={uj-\y+), 



!..+ ) = |^-)F, (c.-HF(c.+ |. 



These spectral data have the properties 

F^=F, [io,F]=0, 



0. 



(2.34) 
(2.35) 

(2.36) 



In terms of the (p, q) representation, the last equality means that there exists a function 
/ of two real variables such that we have the representation 



F(p;q) = 5{p2 - Pi{pi + 2qi^))f{pi +qm;qisR), 



(2.37) 



and then ( 2.36 ) means that f{a ,f3) — f{j3,a). In these terms, we can use ( 2.21 ) to 
rewrite the first equation in ( 2.35 ), for instance, in the form 



X+(x,A:) = / dax"(a;,a)e"(^('=)-^("»/(a,A:). 



(2.38) 



We thus see that Eqs. ( 2.35| ) are equivalent to the standard (see |1^, Q) formulation 
of the nonlocal Riemann-Hilbert problem for the Jost solutions of the NS equation. In 
what follows, we assume the unique solvability of these equations. 



3 Integrable symmetries and time evolutions 

3.1 Deformations of the spectral data 

We now consider two potentials: u{x) with the corresponding operators L, M, |z^), and 
(w| and the spectral data F and the potential u{x) with the corresponding L, M, |j>), 



{uj\, and F. Then the two resolvents are related via (2.28) and (2.30) using "scattering 
on a nontrivial background" pOl , 



where we introduce an operator rj such that 

ri=\V){uj\, r;"^ = |j/)(w| 



This operator is unitary; because of ( p.27[) 



v^ = v ^; 



(3.1) 
(3.2) 
(3.3) 



its kernel in the (p, q) representation does not depend on q2, r]{p;q) = ri{p;qi); it is 
analytic for qig ^ and satisfies the asymptotic condition 



lini r7(p;qi) 

qi^oc 



S{p). 



(3.4) 



For t he boundary values of 77 on the real qi axis (defined by analogy with ( ^.31 ) and 
( ^.32| )), we obtain 

,7± = |?±)(c.±| (3.5) 



from (3^), and from ( p.35| ), we then find 

?7+-?7" = \i?-){F ~ F){uj+\ 



(3.6) 



We suppose that a curve in the spectral-data space is parametrized by t. Corre- 
spondingly, all objects M, L, U, |i^), and (w| become dependent on t; we therefore have 
a symmetry of the NS equation. Let the nontilde variables correspond to some value t 
and the tilde variables correspond to the value t + t. Then 77 — ri(t,T), and by ( ^.28 ) 
and (U), 7?(i,0) =/. We set 






(3.7) 



T=0 



This operator is Hcrmitian by ( |3.3| ), A'^ — A, and its kernel v4(p, q) is analytic in the 
upper and lower half-planes and by ( |3.4[ ) satisfies the asymptotic condition 

lim ^(p;qi) =0. (3.8) 

It follows from (p]q) that the boundary values of this operator at the real axis satisfy 



f)F 

A+ -A- = -i\v-)^{uj^ 



(3.9) 



It is clear that the analyticity properties of A(p;q), together with (3.8) and (| 
uniquely determine the operator A in the (p, q) representation via the Cauchy formula 



A(p;q)^- 



qiK - qi V 9t 



(3.10) 



we recall that the kernels of all objects in the r.h.s., \u), F, and {uj\, depend neither on 
q2 nor on qijj in the (p, q) representation. We can use ( |2.4| ) to rewrite this kernel in the 
X representation as 

Aix,x';q) = -zsgnq,e-^^(-^-<^eiix,-x[)q,) ('|^-)^(c.+ |^ (x, x'; g), (3.11) 

where the last factor is independent of q by construction. 

If A is known, then for the resolvent and the operator L itself, we have 



dM 



i[A,M], 



dL 
'dt 



i[A,L] 



(3.12) 



from ( |3.l[ ) and (pj). Moreover, using (2.28) to rewrite (3.2) as \v) = 'ij\v) and {uj\ — 
{uj\rj^'^ , we obtain the derivatives of the Jost solutions along the curve from (pj): 



iA\v), 



-i{uj\A. 



(3.13) 



dt ' " dt 

The flows of the operator L, the resolvent, and the Jost solutions are determined 

using the 1 ST a nd a re gi ven by the operator A, i.e., by dressing the flow of the spectral 

data with ( 3.10 ) or ( 3.11 ). Taking into account that Lq is independent oft by definition, 

we obtain 

dU 
z— = [Ai] (3.14) 

for the potential U from ( |2.8| ) and ( 3.12| ). Due to ( ^.10 ), the kernel of the l.h. s. is 



independent of q in the (p, q) representation, and A{p\ q) decays as qi — > oo. Taking (2. 



and (2.9) into account, we see that commutators with zl?2 and U in the r.h.s. also decay; 
the nonzero result therefore gives only the commutator with D\, 



3JJ 
^—^-^l^A-2{^,^A)D,, 



where we introduce 



d:A = [D, 



[D,.A] 



(3.15) 



such that {dxjA)(x,x' ,q) — Axj{x,x' ,q) + Ax'Xx,x\q) in the x representation. Again, 
d^ A decays as qi ^ cxo; therefore 



— -29 A(-i) 
dt ~ ^' 



(3.16) 



where A^ ^^ is the operator whose kernel A ^{p, q) is equal to the residue of A{p; q) at 
infinity, i.e., 

^(-i)(p;q)= lim A(p;q)qi. (3.17) 



From (3.10), we have 



A(-l)(p;q) = i-|dq;^(^|^-)^(,.+ |^(p;q') 



(3.18) 



and from ( 3.1l| ), 



Aix,x';q) ^ S{x, ~ x[) (\,^-)^{u;+\\{x,x'-q) (3.19) 

in the x representation. By construction, the last factor is proportional to d{x2 —x'2) and 
is independent of q. We can therefore see that the r.h.s. of ( |3.16) is i ndeed a multiplication 
operator in the x representation, as it must be according to (2.10). 



3.2 Linear symmetries of the spectral data 



To obtain integrable evolution equations associated with the NS equation from this con- 
struction, we must consider Hnear symmetries of the spectral data F. We therefore 
impose the condition that there exist operators a such that 



dF 



a^F ~ Fa^ 



(3.20) 



Because of (2.37), we can chose kernels of these operators in the (p, q) representation 
that are independent of q2 and qiQ, i.e., without loss of generality. 



Properties ( ^.36 ) of the spectral data show that we can impose the conditions 

again without loss of generality. 

Returning to (3_^) and using ( 2.35 ) and ( 3.20| ), we obtain 



A+ - \j^+)a+{uj+\ =A-- \iy-)a-{uj-\. 



(3.21) 



(3.22) 



(3.23) 



By construction, A^ , 



and 



have analytic continuations in the upper and lower 



half-planes. It is therefore natural to impose the condition that there exists a function 
a(p;q) (analytic or meromorphic) in qi for qig ^ and independent of q2 such that 



a=^(p;q) 



lim a(p;q). 



(3.24) 



We can thus consider different symmetries described by operators a with kernels a(p; q) 
that have singularities of different types in the complex plane. Here, we mainly consider 
the case where a(p; q) depends on qi polynomially. Because, as is well known, the 
expansion coefficients at infinity of the Jost solutions of the NS equation are given by 
recursion relations in terms of the potential u, we call such symmetries explicit. In 
particular, because a has no discontinuity on the real axis for such symmetries, i.e., 



conditions ( 3.22| ) give 



[L„a]=Q. 



We then have 



dF 



instead of ( 3.20| ), and ( 3.1Cl|) can be written as 



A^{W)a{u;\)^, 
P ^ {W)a{u;\)^, 
\iy)a(Lu\ ^A + P, 



(3.25) 
(3.26) 
(3.27) 



(3.28) 
(3.29) 
(3.30) 



where (•)+ and (•)_ denote positive (polynomial) and negative parts of the 1/qi expan- 
sion at infinity of the kernels in the (p, q) representation. From relations ( 3.28| )-( 3.30 ) 



with (2.28) and the commutativity of a with Lq and with Mq taken into account, we 
obtain 



.^^IP.Ml, 



dt 



^^ = iP^L], 



P\u) - \i^)a, 



dt 

d{Lj\ 



a{uj\ — {uj\P 



(3.31) 
(3.32) 



instead of ( 3.12| ) and (3.13). 

It is easy to prove that all polynomials a with properties ( 3.26D can be obtained as 
linear combinations of the operators 



«™,„ = ^{^r,A"}, m,n>0, m + n>0. 



(3.33) 



Here, 



A = Xi +2^X2151, 



(3.34) 



and we introduce the operators Xj that are operators of multiplication by Xj in the x 
representation, 



XJx,x';q)—XjS{x — x'), j = 1,2. 



It is easy to verify that 



X] = X^, [D^,X,]^5^^„ 

At=A, [D,,A]=I, [D^,A]=2iD„ 

[Lo,A]-0. 

Symmetries generated by the operators am,n were called additional in 
consider time evolutions associated with their subclasses. 



(3.35) 

(3.36) 
(3.37) 
(3.38) 

We here 



3.3 Subsets of commuting symmetries and hierarchies of evolu- 
tions 



In general, the flows introduced in ( 3.33| ) do not commute. On the other hand, it is clear 
that in the linear span of these operators, there exist subsets of mutually commuting 
operators. We consider two operators a and a', let t and t' be the corresponding evolution 
parameters of these symmetries, and let A, A' and P, P' be defined by a and a' as in 
( [3.2g ) and ( 3.29| ). In terms of the spectral data, the commutativity condition for these 
flows is simply 

[a,a']=0, (3.39) 

where we assume that a and a' are independent of t and t' . In terms of A and P, this 
condition is more involved. Indeed, from (^.32), 



d_ 

at 



WWH 



A,W)a'{cu\ 



and vice versa; we therefore have 



^ {\,y)a'{Lu\) - ^{W)aH) = ^W)[a,a']{uj\ - t[A\A]+z[P',P]. 



dt 



(3.40) 



10 



In case (3.39), the polynomial and the decreasing parts of this equality can be separated, 
and we obtain the Zakharov-Shabat conditions 



dA 



dA' 
~dt 



z[A,A% 



dP 



dP' 



l[p,p'] 



(3.41) 



As we see from (3.33|) and ([3.34), am,n is a differential operator of order m + n. If 



TO and n are not relatively prime, i.e., if there exist numbers fc ^ 1, 



and n' such 



that m = m'k and n — n'k, then am,n is equal to {am\n')'' pl us a l inear combination of 
some lower «„//_„", to" + n" < m + n. The linear span of set ( 3.33| ) is thus decomposed 
into the union of subsets of commuting flows. These subsets are labeled by pairs of 
relatively prime numbers {k,l), the operator of the lowest order in the {k,l) subset is 
akj, and all elements of this subset are equal to {ak,i)"^, to, = 1, 2, . . . (or to their linear 
combinations). Because these subsets of mutually commuting operators are disjoint, we 
can consider only symmetries belonging to one (fc, I) subset. For simplicity, let a denote 
0'k,i- We thus study symmetries determined by the operators a™, to. = 1,2,..., and let 
the corresponding evolution parameters be denoted by t^, 



^^^[a"\F], to = 1,2,, 



(3.42) 



To emphasize the commutativity of these symmetries, we use the term times for these 
parameters. 



Using ( 3.42 ) and inverse problem ( 2.35 ) for the Jost soluti ons, we ca n int roduce 
dependence on these times for all objects of the theory; thus, by ( 3.31 ) and ( 3.32| ), 



dM , 



■ dL 

.'9(^1 m/. .1 



dt. 



(lo\P„ 



whereby (|l29|) . 



P„ = (k)a™(^|)_ 



(3.43) 
(3.44) 

(3.45) 



To be more precise, we mark the original operators with a hat: F, U, jP), (a)|, etc. Then 
the t-dependent operators F{ti,t2, ■ ■ .), etc., are determined by the corresponding initial 
data, 

^(ii,i2,---)lti=*.=„,=o = i^, (3.46) 

etc. We emphasize that all time variables are here introduced in addition to the original 
variables of the kernels of all operators participating in the construction of the spectral 
theory of the NS equation. It is clear that the time variables are distinct from the 
variables of the kernels of operators; in particular, the derivatives with respect to these 
sets of variables commute. 

The times tm are naturally ordered by powers of the differential operators a™. We 



mention that evolutions ( 3.42) can also be considered if a{p; q) is a meromorphic function 
in the complex domain of qi . In this case, of course, Pm stands for the singular part of 
the Laurent expansion of |i^)a'"((jj|; it is then given in terms of not the potential U but 
the values of the Jost solutions at the poles of a. It is important that in all these cases 
independently of the exact choice of operator a, evolutions ( ^.42| ) are associated with the 
KP equation and its hierarchy. Indeed, for arbitrary operators a and F, there exists the 
simple commutator identity 



Ma,F]]~-[a^ 



[a^F]] + -[a, 



.FM^O, 



(3.47) 
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which imphes that the spectral data F satisfy the differential equation 

d^F 3 d'^F 1 d^F 



dt^dt^ 



Adtl 



4 dt\ 



= 0, 



(3.48) 



i.e., the linearized KPI eq uatio n, and so on for higher linearized equations. 
Moreover, because of ( 3.44 ), we have the operator equations 






PiW) 



\v)a, 



(3.49) 



d'^\v) ( o dF\ , , dlv) 






P. 






P? 



w + 



,. d{PiW)) 



which, in the case where dPi/dti ~ 0, are simplified to 



.d\i 



dU 



.d\v) . 

Otr, 



dtl 

d'^\v) 

dtl 



= {P2 - Pi) 1^) + 2* 



d\i 



dt^ 



(P, ~ P!) \i^) + iiP, 



9H 



(3.50) 
(3.51) 



here, the first equation is the obvious analogue of ( 2.25 ) with the operat ors P ^ — P^ and a 
playing the role s of t he potential and of the spectral parameter. Using ( 3.4g| ) and ( 3.50|) , 
we can rewrite (|3.51 ) as 



d\v) ii d^\v) ld^\v) 



dts 2 dtj^dt2 



2 dtl \ ^ 2 ^ 2 ^ ^ 



l'^) 



2 ^ 



.d\u) d^\u) 



dt^ 



dtl 



The derivative of this equation with respect to ii after -^a is replaced using ( |3.50D gives 
the "dressing" of Eq. $M) 



■\v) id''\v) ^ \d'^\v) 
^^4 dt\ 



dt^dt^ 



4 dtl 



Pz + \p!~\{Pi^P2}.Pi 



iidP^ 



- M ^3 + \p! - \P2P1 



i(^^ 



Pi 



d\y) 



d\v) 
'di~ 



dt^ dtl 



(3.52) 



where we use ( |3.4lD and {•, •} stands for the anticommutator of operators. 

The simplest examples of polynomial fiows of type (3.42) are given by two special 
subsets (1,0) and (0,1), which, by (3.33), are respectively 



aifl = iD^, 
ao.i = A. 

We consider these evolutions in more detail below. 



(3.53) 
(3.54) 
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4 Examples of the simplest evolutions 

4.1 Times t„n 



We first briefly demonstrate the above technique applied to the standard flows (3.53). 
The corresponding times and polynomials are denoted by t„j g and Pm,Q- We thus have 



dF 



dt 



m,0 



■ d\v) 
dt 



[{iD,r,F], m=l,2, 



P^M - \v){^D,r 



m.Q 



P. 



i\^.){^D,y^co\l 



(4.1) 

(4.2) 
(4.3) 



Using representation ( 2.37 ) for the kernel of operator F, we obtain time evolutions 
of the spectral data in the standard form 



a/(a,/3) 



dt 



m,0 



zKa'" -/?"]/(«, /3). 



(4.4) 



To calculate the polynomials Pm,o explicitly, we use the asymptotic expansions of the 
Jost solutions 

OO C30 

\i,) = J2 '^n{2iDi)-'-, {u;\ = ^(2iDi)-"o;„, (4.5) 

n=0 n=0 

which have the exact meaning in our approach of the 1/qi expansion at infinity for 
|j/)(p;q) and the l/(qi + Pi) expansion for {uj\{p;q). It is known that the coefficients 
of these expansions depend on the half-plane of qi, but here we can ignore this fact, 
referring to |0 and |18| for details. To calculate these coefficients, we can use either the 
integral equation for the Jost solutions or bilinear relations (2.28|) and the relation 



U = iD^ + Dl - \j^)(tD2 + Dl){, 



LJ 



which follows from (2 



(p. 91) , and ( p. 29]) . We then obtain 
j/i = -cji, d^^iy^ = ill, 

.2 ,, , 1q .,2 



dx,^2 = ^X2^1 - «5:ri^l + -j^x^V^, 



(4.6) 

(4.7) 
(4.8) 
(4.9) 

(4.10) 



and so on, where notation (3.15) is used. We omit the corresponding recursion relations 
because we do not need the highest coefficients here. Using these formulas and (4.3), we 
obtain the combinations involved in (3.49), (3.5C), and ( 3.52| ): 

(4.11) 
(4.12) 

(4.13) 



^2,0 ~ A,0 — '^' 



For the first expression in (3.52), we have 



p. 



3,0 



1 3 

2 "l,0 ~ Tl "2,0' "l.ol' "l,0 



3i 9P2,o Si dU 



Kv 



(4.14) 
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By (3.49) and (3.5C), \v) thus satisfies the equations 

d\v) 



.d\v} d'^\i') 



dt 



2,0 






"'■1,0 



(4.15) 
(4.16) 



which, due to ( 



U), ( ^^ , (|1|), and dZTsI ), are exactly 

a|zy) 



9^, 



-a,J^). 



(4.17) 



Correspondingly, from (4.5) and (4.8), we obtain the equalities for the potential U 

dU „ _ dU 



dt -^-^' 

"'-1,0 



at. 



a.„[/. 



Using representations (p.lC|) and (2.21), we obtain 



dxjx, k) dxjx, k) 

—qI ^Xa:^{x,k), -^^ =X:.,(a;,A:), 

du{x) 



Ux^ix), 



5*2,0 

du{x) 



UxS^)- 



(4.18) 

(4.19) 
(4.20) 



Bf 

These equalities are often considered the reason to identify ii^o — xi and t2.o — X2- Here, 
as mentioned in discussing (3.45)-( 3.46| ), we write 

X(''1,0' *2,0Fl' ^2; k) = X[Xi + tiQ, X2 + ^2,0' k), 

"(*i,0: t2,o\xi,X2) = u{x^ + ti 0, X2 + t2,o)- (4-21) 

Now, we see that due to ( [4.1^ ), expression ( 4.14 ) is equal to zero; ( 3.52 ) is therefore 



reduced to 



9*1,0^*3,0 " 4 af;; + 4 9^ 



= ■:; (d.^U ^ d^X, 



'9*1,0 



3U .dM 52|j/) 



dt 



2,0 



dti, 



and from (4.5) and ([4.8|), we obtain the standard KPI equation 

d'^U 3 d^U 1 d^U 3 d^U"^ 



5*1,0^*3,0 4 9t2_„ 4 9tf,o 4 9^2^0 



(4.22) 



(4.23) 



4.2 Times t 



0,m 



In this section, we investigate evolutions determined by (3.54). The corresponding times 
and polynomials are denoted by to^m E^nd Po,mi ^iid we again emphasize that derivatives 
with respect to to,m commute with the x variables and the times t^.o considered above 
are not switched on now. Here, we thus have 

■ ^^ -[A",F], m=l,2,..., (4.24) 



dt 



O.rn 



R .|^)-|z.)A^ 



"*0,m 



(4.25) 
(4.26) 
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where the operator A is defined in (3.34). Time evolutions of the standard spectral 
data (2.37) are now 



dfia,P) 



■A 

da 



d_ 

"d(3 



/(«,/3), 



(4.27) 



and for simplicity in what follows, we consider only m — 1,2,3. 

The polynomials Po.m can be explicitly written using ([4.5[), (4.7), and ( [4.10[) ; it must 
be taken into account that X2 commutes with |i/) and {uj\ because of ( p. 19 ), ( 2.20D , 
and (|3.35). For the combinations involved in (3.49), (3.50|), and (3.52|), we thus obtain 



PoA = A, 

Po,2 - Po,i = 4X|C/, 



0,3 -r ^^0,1 - ^^0,2^0,1 = G^Xlid^^U- d,,^v^) - QiXld,,^{X,v,), 



Pr 



and for the first expression in (p.52|), we have 



(4.28) 
(4.29) 

(4.30) 



Pc 



1 3 3 

0,3 T" 2 '''I ^ 7 I 0.2j "0,1/1 "0,1 



= SiX^ 



dU 
dto,2 



3^ 5Po,2 ^ 
4 9io 2 

-4X^[d,^{X,U)+d,^{X,U)]]- (4.31) 



Then, by (|1|) and (|3^ ), 



= -t[A,\,,)] = -t[X„W)] +2X,d,Ji.), 



dliy) d'^\v) „ d\v) 
1 ! L J \ L — On \ L 



A + 4X|C/|i.), 



(4.32) 
(4.33) 



and the second equation, i.e., the NS equation with respect to these times, is thus the ex- 
act a nalo gue of ( 4.16 ) with the role of the potential played by the multiplication operator 
(cf. (S) 



W^^XlU, W(x,x ,q)^w{x)8{x-x), w{x)=4x\u{x). 



Expanding these equations by ( [4.5[ ) and using (4.7) and (4.8), we derive 



^=2XA^U, 

'^''0,1 



at 



0,2 



= AX,[d,^{X,U) + d,^{X^U)], 



(4.34) 



(4.35) 



where we take [Xi, |i^)] ~ D^^ into account 
The evolution with respect to 
1.31) is equal to zero due to 



The evolution with respect to io,3 via (3.52) is again simplified because the expression 
in (4.31) is equal to zero due to (4.35). From (4.29) and (4.30|), we then obtain 



d\u) 



dto,idto,3 4 dtl, 4 dtl. 



%Xl {X.id^^U- d^^u,) - d,,^{X,y,)) ^ + 



''0,2 



^''0,1 



(4.36) 
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Using the above relations, we can extract the term of order D^ ^. For the potential W 
in Eq. ( 4.33D , we then obtain exactly the KPI equation in terms of the times io,i, ^0,2, 
and to, 3: 

d'^W 3d^W 1 d^W _ 3 d^W^ 

'^'■0,l*^''0,3 ^ "''0,2 ^ '^'^0,1 ^ '^''0,1 

4.3 Transformation of dependent and independent variables 

As shown in the previous section, taking the time io.m into account gives rise to a new 
dependent variable w{x) via (4.34). Because of ( 4.37| ), this function of five variables, 
^^(^0,11 io,2, ^osl^;), satisfies the KPI equation with respect to the io,m variables. 



d'^w 3 d'^w 1 d^w 3 d'^w'^ 

dto,idto,3 ~ 4 9^ + 4 <; = 4 9^' 

for all X, and Eq. ( [4.35D , together with l \2.1(\ ) and ( [4.34| ), implies 



(4.38) 



— = ^x^w^^ (a;), —^ = Ax2[x^w{x)^^ + x^wix)^^ 



dU 



(4.39) 



which suggests a transformation of the independent variables as well. Indeed, let 

-1 



Xl 



4a;, 



(4.40) 



then instead of ( 4.39 ), we obtain 
dw{x) 



dw{x) 
Wz^[x), — ^w,^(x). 



dt 



(4.41) 



0,2 



This means that to,i and io,2 shift the variables zi and Z2. In other words, if we introduce 



'P(*0,l'*0,2k) ='^(Vl'V2l^)' 

then (cf. (p^)) 

'^(Vl'*0,2|2) = ^(^1+ Vl'^2+ V2); 

letting w denote the value of w at io.i = ^0.2 = 0, we obtain 



(4.42) 



w{x) = ip 



Xi — 1 



and 



"' ' \ I-4X2V1 1-42:2*0.1 



(4.43) 



(4.44) 



We see that the class of initial data and, correspondingly, the class of solutions of 
the KPI equation can be essentially extended in this way. Indeed, we seek a solution 



■y^Cio,!! ^0,2,^0,3) of Eq. (4.38) that satisfies the initial data 

ipito,l,to,2,0) = <?(io,i,io,2)- 
Let 

■^(^0,lJ*0,2i^0,3|2:) = 95(^0,1 + ^li^0,2 + ^21^0,3)7 



(4.45) 
(4.46) 
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where z is given by ( }4.40| ). This function also solves ( |4.38| ), and by ( |4.45| ), Eq. ( [4.43[ ) 



pves w{x) — u;(0, 0, 0|a;). Switching on the times io,m and constructing the function 
M(io,ij*o,2,io,3|2:) with the above procedure, we conclude from (4.34) and ( 4.43| ) that it 
satisfies 

Transforming the function w(to,i,io,2, ^o.sl^;) — 4a:|u(io,i,io,2, ^o.sl^;) with the transfor- 



mation inverse to (4.40) and setting z = 0, we obtain the solution of the Cauchy problem 



for ((5(to,ii ^0,2: ^0,3) from (4.45). By (4.47), the equation is thus solved for the initial data 



^(*0.1'^0,2) = 7T2-" 



4/2 10^ '4+ / ' 

where u as function of its arguments — not (p — must satisfy the smoothness condition, the 
decay-at-spatial-infinity condition, and the small norm assumption in order to guaranty 
the applicability of the 1ST. 
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